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Fluctuations of the thermal or classical component of the van der Waals force between two dielec- 
tric slabs, modelled as an ensemble of polarizable dipoles which interact via the usual electrostatic 
dipole-dipole interaction, are evaluated. In the model the instantaneous force is a deterministic func- 
tion of the dipole configurations in the slabs and its fluctuations are purely due to dipole fluctuations 
(no background thermal fluctuations of the electromagnetic field are considered) . The average of the 
force and its variance are computed. The fluctuations of the force exhibit normal thermodynamic 
scaling in that they are proportional to the area of the two plates, and even more importantly, do 
not depend on any microscopic cut-off in the theory. The average and the variance of the thermal 
van der Waals forces give a unique fingerprint of these fluctuation interactions. 
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INTRODUCTION 

The Casimir force pQ is the force arising between ob- 
jects placed in a quantum and/or thermal field due to the 
modification of the fluctuations of the field by the pres- 
ence of the objects. As the force is due to a fluctuating 
field, the force itself should fluctuate. What is normally 
given as the Casimir (quantum or thermal due to the 
zero frequency M atsubara mode) force is its average value 
measured in thermodynamic equilibrium. Deriving a fin- 
gerprint of the average force and fluctuations around the 
average force could aid identification of the Casimir in- 
teraction component in an otherwise complicated exper- 
imental setup, showing long-range interactions of mixed 
origins [2]. 

The first analysis of Casimir force fluctuations pertains 
to the quantum context [3j[4]. In this case, despite the 
fact that the average value of the force is finite, the vari- 
ance of the force exhibits an ultra-violet divergence that 
is eliminated using the experimental fact, that the force 
is always averaged over a time corresponding to the tem- 
poral sensibility of the experimental apparatus. From 
this analysis, while the average force on a single mirror 
is zero, its variance is non-zero, due to differences in the 
electric field on either side of the mirror. These fluc- 
tuations of the Casimir force can in turn be related to 
the force exerted on a non-uniformly accelerating mir- 
ror [5]. In addition, the fluctuations of the radiation 
pressure exerted by a laser beam on a conducting sur- 
face can also be analyzed [6 , in this case calculations 
using the electromagnetic stress tensor can be confirmed 
by a kinetic like approach based on photon number fluc- 
tuations. Casimir force fluctuations between perfectly 
conducting mirrors at finite temperature have also been 
examined [7] and it was shown that the high temperature 
force could be derived from the classical Rayleigh- Jeans 



distribution, while the force fluctuations required a full 
quantum treatment before taking the classical limit. In 
all the above cases average forces and their fluctuations 
are seen to arise via the boundary conditions imposed 
by the conductor on the electromagnetic field. In [8, 9 
fluctuations of the Casimir-Polder interaction between a 
polarizable atom and a perfect conductor were analyzed, 
presenting a conceptual departure from previous studies 
as force fluctuations can be related to a physical property 
of the atom, its polarizability, going beyond descriptions 
of objects in terms of ideal boundary conditions. 

In [10] the fluctuations of the thermal Casimir force 
due to a free mass-less scalar field theory with Dirich- 
let boundary conditions on parallel plates was consid- 
ered. The leading term in the variance for two plates of 
area A separated by a distance L was shown to posses 
a cut-off dependent limit (f 2 ) c ~ (/c#T) 2 A/a 4 , where 
T is the temperature of the system, ks is Boltzmann's 
constant and a is a microscopic (ultra-violet) cut off. 
The average value of the force in this system behaves 
as (/) ~ ksTA/L 3 . Thus in terms of the extensive vari- 
able A, the fluctuations of the force A/ obey the usual 
thermodynamic scaling A/ ~ y/A, compatible with the 
notion that local fluctuations of the force at distant re- 
gions of the plates are uncorrelated. Similar results have 
also been found for the fluctuation induced forces on in- 
clusions, such as proteins, in membranes [TTj . 

More recently [12 the fluctuation of the Casimir force 
for scalar fields was examined in a parallel plate piston 
cylinder geometry, with the result that (f 2 ) ~ 2(/) 2 , 
which is clearly at variance with the thermodynamic scal- 
ing found in [3, 4, 10 . Whether thermodynamic scaling 
should hold is far from obvious as all the fluctuation in- 
duced interactions in the studies above are long range 
(corresponding to mass- less field theories). 

Connected with the question whether fluctuations in 
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FIG. 1: (Color online) Schematic drawing of a classical dipole 
model of two dielectric half- spaces at separation L. The slit 
region is a vacuum. Fluctuating thermal dipole configurations 
(shown bold and faint) create a fluctuating thermal Casimir 
interactions that for large L decay to zero strictly with a prob- 
ability one. 



Casimir forces scale thermodynamically, is the intriguing 
existence of fluctuating forces on isolated bodies, depend- 
ing on a microscopic cut-off, as such forces should induce 
movement /diffusion of small inclusions and perhaps also 
fluctuation-induced drag forces, which have been pre- 
dicted in a number of quantum [13 j and thermal [14] 
situations. The theoretical results described above point 
to the existence of regimes, notably in microscopic sys- 
tems, where fluctuations of measured Casimir forces of 
both quantum [15] as well as thermal or critical Casimir 
provenience [16], can become large and thus may be ex- 
perimentally measurable. 

In [17] a model dielectric was introduced, based on a 
continuous polarizable dipole field, showing thermal fluc- 
tuation forces identical to the zero frequency Matsubara 
mode van der Waals forces of the standard Lifshitz the- 
ory with media of the same dielectric constants. Within 
this model it is possible to analyze straightforwardly how 
the van der Waals interaction arises from the correlations 
between dipoles in opposing slabs and how the van der 
Waals force evolves temporally when switched on from 
zero, that is for initially uncorrelated slabs, to its final 
equilibrium value. 

Here we use the very same model to study the fluctua- 
tions of the dipole-dipole induced thermal Casimir force 
in the problem. Because we analyze direct interaction 
between dipoles, it is clear that as the distance between 
the two slabs is taken to infinity, the force itself must go 
to zero. In our model we only take into account the elec- 
tromagnetic field generated by the dipoles in the slabs 
and ignore fluctuations of the electromagnetic field due 
to thermal photons (which cannot be properly included 
in the classical limit considered here [7]). In this situa- 
tion therefore, we do not expect a bulk cut-off dependent 
force fluctuation of the type mentioned above. We find 
that the average force scales as (/) ~ AksT/L 3 while 
the variance behaves as (f 2 ) c ~ A(kBT) 2 / ' L A and tends 
to zero for infinitely distant bodies. We note that the 



variance obtained behaves like the L-dependent part of 
the variance found in [TO] . 

POLARIZABLE FIELD MODEL 

We consider a classical model of interacting dipoles in- 
troduced in [17] . Here we have two slabs of material sep- 
arated by distance L in the direction e z and the Hamil- 
tonian for the system is given by 
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where p(x) is a local dipole field. The interaction is given 

by 



Aij (x, x : L ) — 



(%£(x-x') 
X(x) 



A,(x,x';L), (2) 



where x(x) is the local polarizability of the dipole field at 
the point x and Dij is the usual dipole dipole interaction. 
In slab S{ we set x — X% (f° r & = 1 and 2) and we use a 
coordinate system such that the points in S\ are in the 
half space V ~ : z < and the points in the slab S2 are in 
the half space V + : z > 0. The dipole-dipole interaction 
is then given by 



^° ) (x-x , ) = V l V'G (x-x') 



(3) 



where V indicates the gradient with respect the the co- 
ordinate x, V' the gradient with respect to x' and Go is 
the vacuum Green's function obeying. 



e V 2 G (x-x / ) = -^(x-x / ). 



(4) 



With this notation the dipole-dipole interaction is 
given by 



Ai(x,x';L) = < 
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Note that we have chosen the slabs to be separated 
by vacuum, so that there is no bulk pressure associated 
with the intervening dielectric medium, independent of 
the plate separation L. The thermal Casimir force usu- 
ally given for such systems equals the plate area multi- 
plied by a disjoining pressure, i.e. the difference between 
the confined and bulk pressures, which tends asymptoti- 
cally to zero as the plate separation increases. 

For a fixed configuration of dipoles (where both their 
relative position and orientation are fixed) the instanta- 
neous force on S2 is given by 
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as only the interaction between dipoles in different slabs 
depends on L. Of course each dipole feels both the elec- 
tric fields from dipoles within the same slab and those 
in the opposing slab. However, when we infinitesimally 
displace S2, the whole slab moves and each dipole is dis- 
placed by the same amount leaving their relative sepa- 
ration and consequently their energy of interaction the 
same. Let Si denote the half space z < and S2 the half 
space z > L, so that the force is given by 



/ 



/ dxdx' 2^ Pi ( x ) 

/SixS 2 



_d_ 

dz' 



D$\*-*) 



Pj (x'). (7) 



This expression for the force has an obvious physical in- 
terpretation, which could in fact have been used as a 
starting point for its definition. The electric field due to 
the dipoles in Si at the point x' in S2 is given by 
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(8) 



and the force in the z' direction on a dipole at x' in the 
S2 due to the dipoles in Si is thus 



F(x') = -5>(x')^£y(x') 



(9) 



so Eq. (p| is just f = J s dx. f F(x!). Now writing the 
dipole-dipole interaction in terms of the free Green's 
function we find 



/ = - / dxdx ; V>(x) [ A y^-G^x - x') 

^ixs 2 *-f yoz j 



(10) 
It is important to emphasize here that in this model the 
electric field is a fixed function of the dipole configura- 
tions and that we do not consider additional thermal 
fluctuations of the electromagnetic field. 



FORCE FLUCTUATIONS 



In principle we can compute moments of the force from 



Eq. (10), but the calculation can be considerably simpli- 



fied by considering 
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(11) 
where, as in Eq. (|7|), we have simply replaced the deriva- 
tives in L by derivatives in z' . In order to compute 
the partition function for this system and the correla- 
tion function for the dipole field, we first introduce the 
generating function 



Z[u} = 



TTd[pi(x)]exp -/3H + dx V"u;(x)p;(x) , 
\ JStUSi ■ / 



(12) 



which can be rewritten introducing an auxiliary field <p 
that decouples the dipole-dipole interactions. We can 
now integrate over the dipole field to obtain 
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(13) 



where by setting <\> = —i^), ip can be identified as the 
electrostatic potential. Integrating over the dipole field 
yields the partition function, while the no source term is 
immediately recognizable as the zero frequency Matsub- 
ara mode or thermal (van der Waals) Casimir interaction 
of the standard Lifshitz theory [19] 



Z[0] = I ^(x)]exp (-£ J dx e(x)[V0] ; 



(14) 



with dielectric constants e(x) = eo +Xi when x is in Si or 
S2 and e(x) = eo for x is between the slabs. Written in 



the form of Eq. ( 14 ) , but also from the basic model above, 



it is clear that our model is valid for high temperature 
dielectric systems, retardation effects are neglected and 
also no effects due to conduction electrons is present. 

In the absence of sources, the average value of the 
dipole field is (pi(x)) = 0, and its correlation function 
at non-coinciding points is 



(Pi( x )ft( x 0) 



x(x) X (x / )V,V , ? G(x,x / ), (15) 



where G is the slab geometry Green's function obeying 
V • e(x)VG(x, x 7 ) = -J(x - x 7 ). (16) 

The average force is thermodynamically given via 
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To compute the variance of the force one notices that 
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which can be rearranged to give the Gibbs lemma [18] 
type result 
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The first term on the right-hand-side is easy to calculate. 



The second term is obtained from Eqs. (11) and (15) as 
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(20) 



This can be further simplified using the Fourier-Bessel 
transform of the Green's functions in the r = (x, y) plane 
to get the remarkably simple formula 
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The lateral Fourier transform of the free Green's function 
Go is given by 



G (k, z-z f ) = -^- eM~k\z - z'\) (22) 

2eok 



while 
G(k,0,L) 



2eoexp(— kL) 
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(23) 
where A^ = (e^ — eo)/(Q + eo). This finally yields 

where Li s is the polylogarithmic function defined by 



Li 8 (z) = £ 



(25) 



The standard result for the average value of the Casimir 
force for this system is given by [19] 
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Putting these results together in Eq. ([19]) finally yields 
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(27) 



where g(z) = Lis(z) — \A±{z). The function g(z) is posi- 
tive for z G [0, 1], so that indeed the variance of the force 
fluctuations is everywhere positive. We see that the scal- 
ing of the fluctuations agrees with the L dependent part 
in [10] , but contains no bulk term dependent on an ultra- 
violet cut-off. 

The ratio of the variance to the squared average force 
is then given by 



*-$-*¥«> 



(28) 



where h(z) = g(z)/Li 3 (z) and h(z) attains its maximum 
value 0.08286 at z = 1. Note the results are symmetric 
under the interchange of the labels 1 and 2 and thus the 
average force on each slab is equal and opposite and the 
force fluctuations on both slabs are the same. This has 
to be the case, as from the definition of the unaveraged 
force it is clear that the forces exerted on each of the 
slabs are symmetric. 



CONCLUSIONS 

We analyzed a simple, classical dipole model, that 
leads to thermal van der Waals interactions with the av- 
erage force predicted by the standard Lifshitz theory. In 
this classical dipole model the force at large separations 
not only goes to zero on the average, but much more 
stringently with probability one. The ratio of the vari- 



ance to the squared average force, Eq. (28), furthermore 



gives a fingerprint of the classical thermal van de Waals 
interaction, differentiating it from possible parasitic ef- 
fects and stray interactions, plaguing experiments in re- 
alistic systems [2]. 

In our dipole model there is no divergent cut-off depen- 
dent bulk contribution to the average force. The fact that 
no other length scale appears means that the functional 



form of the force variance in Eq. (27) can be predicted 



solely on the grounds of dimensional analysis and by as- 
suming the thermodynamic scaling, i.e. A/ ~ \^A for 
the fluctuations as a function of the area of the plates. 
The so obtained thermodynamic scaling is in agreement 
with a number of studies of fluctuations of Casimir forces 
in both thermal and quantum systems [3] [4] [TUJ . The re- 
sults of [10] are for thermally fluctuating fields, as here, 
but find a microscopic cut-off depend force fluctuations at 
infinite plate separation. Nevertheless the average force 
has the same behavior in both models. This difference 
comes naturally from an underlying microscopic defini- 
tion of the model. In our model the electrostatic fields 
are generated by the dipoles themselves, being solutions 
of the Poisson equation for the given dipole distribution. 
The fields so generated cannot generate self-forces on iso- 
lated objects, hence as the two slabs are separated the 
total force on each slab decays to zero and by consequence 
so does its mean and variance. In the model of [10] (and 
indeed its quantum counterparts [3, 4 ), the fluctuating 
field exists throughout space and objects immersed in 
the field are not its sources. The objects immersed in the 
field however modify its fluctuations by the imposition of 
boundary conditions. Crucially however, the fluctuating 
field exists in the presence of a single object and hence 
can exert a force, albeit of mean zero, on the object. 
Intriguingly, the nature of force fluctuations may there- 
fore be relevant to the long standing debate (see [20] for 
a recent discussion) as to the physical interpretation of 
the Casimir effect in terms of a shift in zero point energy 
(field effect) or as a van der Waals effect (due to sources). 
It should be noted that fluctuations of the force acting 
on a single particle is compatible with the drag predicted 
on isolated moving objects coupled to a thermalized fluc- 
tuating field |T31 E]. 

Normally in electrostatics forces are evaluated by using 
the stress tensor [21 that follows from equations of the 
form of Eq. (|9|, in order to compute forces generated 
by electric fields on charges and/or dipoles. The usual 



stress tensor given for dielectric systems depends on the 
local dielectric constant [21] which is already an object 
derived from thermodynamic averaging, relating the av- 
erage local polarization field to the local electric field. 
This means that the validity of the use of the preaver- 
aged stress tensor to compute fluctuations of forces in 
these types of models is not obvious. The stress tensor 
before averaging is a fundamental object derived from 
the force exerted on charges by the electric field. It is 
indeed this stress tensor that should be used to com- 
pute force fluctuations. However, in many coarse grained 
theories used to describe fluctuating systems, dynamical 
variables, which are equivalent to charges or dipoles in 
electrostatic language, have been integrated out in for- 
mulating the theory and it is thus not always obvious 
that they contain sufficient physical information to yield 
correct predictions for force fluctuations. 
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